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EXTREMAL FUNCTIONS IN DE BRANGES 
AND EUCLIDEAN SPACES II 

EMANUEL CARNEIRO AND FRIEDRICH LITTMANN 


Abstract. This paper presents the Gaussian subordination framework to generate optimal one-sided ap¬ 
proximations to multidimensional real-valued functions by functions of prescribed exponential type. Such 
extremal problems date back to the works of Beurling and Selberg and provide a variety of applications in 
analysis and analytic number theory. Here we majorize and minorize (on M. N ) the Gaussian x i— ► e~ 
where A > 0 is a free parameter, by functions with distributional Fourier transforms supported on Euclidean 
balls, optimizing weighted L 1 -errors. By integrating the parameter A against suitable measures, we solve the 
analogous problem for a wide class of radial functions. Applications to inequalities and periodic analogues 
are discussed. The constructions presented here rely on the theory of de Branges spaces of entire functions 
and on new interpolations tools derived from the theory of Laplace transforms of Laguerre-Polya functions. 


1. Introduction 


In this paper we continue our study of the Beurling-Selberg extremal problem and its applications, in 
an unprecedented general setting. We develop here the Gaussian subordination method to produce optimal 
(in weighted L 1 -metrics) one-sided approximations of multidimensional real-valued functions by functions 
of exponential type. In particular, this plainly becomes the most powerful framework to deal with this 
sort of extremal problems, extending a number of previous works in the literature (e.g. mmm in the 
one-dimensional case, and mi in the multidimensional case). 

We follow closely the notation of our precursor m to facilitate some of the references. 


1.1. The extremal problem. Throughout the text we use the standard Fourier analysis notation as in 
[44] . We denote vectors in 1^ or C N with bold font (e.g. x, y, z) and numbers in K or € with regular font 
(e.g. x , y, z ). For z = (zi, Z2, • • ■, zn) G C n we let | ■ | denote the usual norm |z| = (|zi| 2 + • ■ • + [zn] 2 ) 1 ^ 2 , 
and define a second norm || • || by 


||z|| = sup 


N 

^ ^ Zn t n 

n=l 


t e 


„N 


and It I < 1 


If F : C N —> C is an entire function of N complex variables, which is not identically zero, we say that F has 
exponential type if 

t(F) := limsup ||z|| _1 log |F(z)| < oo. 

IM|->oo 

In this case, the number t(F) is called the exponential type of F. Naturally, when N = 1, this is the classical 
definition of exponential type. When N > 2, our definition is a particular case of a more general definition 
of exponential type with respect to a compact, convex and symmetric set K C (cf. (44j pp. 111-112]). 
In our case, this convex set K is the unit Euclidean ball. 
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We address here the following extremal problem: given a function T : 1 N — > R. and real parameters S > 0 
and v > —1, we seek entire functions £ : C N —>■ C and Ai : C N —> C such that 

(i) C and Ai have exponential type at most 6. 

(ii) C and Ai are real-valued on R w and 

£(x) < .F(x) < M(x) (1.1) 

for all x € R^. 

(iii) Subject to (i) and (ii), the values of the integrals 

/ {.Ad(x) — .F(x)} \x.\ 2l/+2 ~ N dx and f {.F(x) — £(x)} |x| 2 " +2_JV dx (1.2) 

Js. N Jr n 

are minimized. 

In the one-dimensional case with respect to Lebesgue measure (i.e. v = —1/2), this problem and its 
applications have a long and rich history, dating back to the works of Beurling and Selberg for the signum 
function and characteristic functions of intervals, respectively (see [SI US]). The one-dimensional theory 
was substantially developed over the past years as we have reached the solution of this extremal problem 
for a wide class of even, odd and truncated functions [H [10j [12l [T3] [23] [33] [34] [33 HU- The fact that these 
extremal majorants and minorants have compactly supported Fourier transforms (Paley-Wiener theorem) 
allows for a variety of interesting applications in number theory and analysis, for instance in connection 
to: large sieve inequalities mmm, Erdos-Turan inequalities Hilbert-type inequalities 

pniii 145] . Tauberian theorems [23] . inequalities in signal processing m , and bounds in the 
theory of the Riemann zeta-function SEuzniainiEDi. Additional interesting works related to this theory 
include 

The weighted one-dimensional extremal problem (including the possibility of other weights rather than 
the power weights) has an intimate connection to the theory of de Branges spaces of entire functions, as 
observed in the remarkable work of Holt and Vaaler m for the signum function. This was later extended in 
[8] for a class of odd and truncated functions, and in QT] for a class of even functions. This extremal problem 
was considered in [5j for characteristic functions of intervals, with applications to the pair correlation of zeros 
of the Riemann zeta-functon (where the relevant weight is {l — ( s ™^ x ) 2 } dx). De Branges’ theory has also 
played an important role in the recent works [29] [35] [36] . 

In the multidimensional case, this problem presents additional difficulties, and has only been addressed 
in two occasions: in the work of Holt and Vaaler 123, for characteristic functions of Euclidean balls, and in 
our precursor m for the exponential function x4e A l x l, where A > 0, and for a family of radial functions 
with exponential subordination given by 

/»00 

x / {e -A l x l — e _A } dcr(A), (1.3) 

Jo 

where a is a suitable nonnegative Borel measure on (0,oo). The class (11.31) includes examples such as 
x4 — log |x| and x 1 —> T(a + l){|x| _a_1 — 1} for —2 < a < 2v + 1. The primary objective of this paper 
is to extend the framework of providing a complete solution of our extremal problem Ot - (HD for 
the Gaussian x 1 —> e —7rA l x l 2 , where A > 0 (that is not included in the class (11.31) 1. and generalizing the 
construction to a wide class of radial functions that strictly contains m- Previously, the special case 
v = —1/2 and N = 1 of the Gaussian subordination framework was obtained in [12] with an interpolation 
method that does not extend to arbitrary v and N. The results here are based on a different approach 
described in Section [2] which is substantially simpler even for v = —1/2 and N = 1. 
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The radial structure is a crucial element in this work, in the sense that it allows us to connect the 
multidimensional problem to a weighted one-dimensional problem, which we solve using the rich theory of 
de Branges spaces of entire functions [5]. 


1.2. De Branges spaces. Before stating our main results, we must recall the basic features of the theory 
of de Branges spaces [3] . Throughout the text we write 

U = {z £ C; Im (z) > 0} 

for the open upper half-plane. An analytic function F : U —> C is said to have bounded type if it can be 
written as a quotient of two bounded and analytic functions in U. Functions of bounded type in U admit a 
Nevanlinna factorization [3] Theorems 9 and 10], from which we draw the number 

v(F ) := limsup y~ l log \F(iy)\ < oo, 

y—y oo 

called the mean type of F. It also follows from [3] Theorems 9 and 10] that the mean type verifies 

v{FG) = v(F)+v(G) 

whenever F and G have bounded type in U. 

If E : C —> C is an entire function, we let E* : C —> C be the entire function given by E*(z) = E{z). We 
say that FI is a Hermite-Biehler function if 

\E*{z)\ < \E(z)\ (1-4) 


for all z £ U, and we say that E is real entire if its restriction to R is real-valued. If E : C —> C is a 
Hermite-Biehler function, we consider the vector space 'H(E) of entire functions F : C —> C such that 

/ OO 

\F( X )f\E(x)\- 2 dx<oo, 

-OO 


and such that F/E and F*/E have both bounded type in U and nonpositive mean type. The space 'H(E) 
is in fact a reproducing kernel Hilbert space , with inner product given by 


{F,G) e 



F(x)G(x)\E(x)\ 2 dx. 


This means that for each w £ C the evaluation functional F <—> F(w) is continuous, and thus there exists a 
function z K(w, z ) £ H{E) such that 


F(w) = (F,K(w,-)) e (1.5) 

for all F £ T-L(E). The function K(w, z) is called the reproducing kernel of T-L(E). 

We denote by A and B the following companion functions 

A{z):= l -{E{z) + E*{z)} and B(z) := l -{E(z) - E*(z)}. (1.6) 

Note that A and B are real entire functions such that E(z) = A{z) — iB(z). The reproducing kernel is then 
given by [3] Theorem 19] 


K(w, z) 


B(z)A{w) - A(z)B(w) 
7 t(z — w) 


E(z)E*(w) - E*(z)E{w) 
2ni(w — z) 


(1.7) 


and when z = w we have 


K(z,z) 


B'(z)A(z) — A'(z)B(z) 
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( 1 . 8 ) 







From m we obtain 


0 < \\K(w,-)\\% = {K(w,-),K(w,-)) e = K(w,w), 

and we notice that K{w,w) = 0 if and only if w £ R and E(w) = 0 (see for instance [27] Lemma 11]). 

1.3. Approximating the Gaussian in de Branges spaces. Throughout the paper we let E be a Hermite- 
Biehler function satisfying the following properties: 

(PI) E has bounded type in U ; 

(P2) E has no real zeros; 

(P3) z i—^ E[iz) is a real entire function (or, equivalently, E*{z) = E(—z) for all z £ C). 

(P4) A,B H(E). 

The fact that E is a Hermite-Biehler function that satisfies (PI) implies that E has exponential type and 
t(E ) = v{E). This is a consequence of Krein’s theorem (see [SD] or EZ1 Lemma 9]). Also, from the Hermite- 
Biehler condition (II. 41) . we note that the companion functions A and B only have real zeros, and moreover 
these zeros are all simple due to (11.81) and property (P2). We are now able to state our first result, the 
solution to a weighted one-dimensional extremal problem for the Gaussian, with respect to a general de 
Branges metric. 


Theorem 1. Let A > 0. Let E be a Hermite-Biehler function satisfying properties (P1)-(P4) above. Assume 
also that 


/ OO 

e-" A| “ |2 \E{x)\~ 2 dx < oo. 

-OO 


The following properties hold: 


(i) If L : C — > C is an entire function of exponential type at most 2 t(E) such that L{x) < e 7rA l x l for 
all x £ R then 

r°° „ _ e-’r-MJI 2 

/ L(x)\E(x )| dx< TqTTT’ ( L9 ) 

J-oo A (£)=0 

where the sum on the right-hand side of Ob is finite. Moreover, there exists an entire function 
z i—^ L(A 2 , A, z) of exponential type at most 2 t(E) such that L(A 2 , A, x) < e -7I ’ A M for all x £ R and 
equality in m holds. 

(ii) If M : C —>■ C is an entire function of exponential type at most 2 t{E) such that M[x) > e -7rA l x l for 
all x £ R then 

r°° „ e -^|c| 2 

/ m {x)\E(x)\ dx> fFTFTf 1 ( L1 °) 

J-CO B(J)=0 

where the sum on the right-hand side of (11.101) is finite. Moreover, there exists an entire function 
z i—^ M(B 2 ,A,z) of exponential type at most 2 t(E) such that M{B 2 ,A,x) > e~f or a ll a: £ R 
and equality in (11.101) holds. 


There are several important ingredients that come together in the proof of Theorem [TJ Among these we 
highlight the following: (a) the construction of suitable entire functions of exponential type that interpolate 
the Gaussian at the zeros of a given Laguerre-Polya function; (b) the decomposition of an L 1 (R, \E(x)\~ 2 dis¬ 
function of exponential type 2 t(E) as a difference of squares of functions in 'H(E)-, (c) the use of Plancherel’s 
identity on the Hilbert space H(E). Since A and B are defined in terms of E , the extremal functions of 
Theorem[l]depend implicitly on E , and the particular choice of notation z i—>• L(A 2 , A, z) and z i-» M(B 2 , A, z) 
will become clear in Section [2] 
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1.4. Approximating the Gaussian in Euclidean spaces. We shall use our Theorem Q] to obtain the 
solution of the multidimensional Beurling-Selberg extremal problem m - m for the Gaussian. As in 
QH123, the main insight here is the specialization of Theorem |T] to a particular family of homogeneous de 
Branges spaces, suitable to treat the power weights \x\ 2lJ+2 ~ N dx. We start by briefly reviewing the basic 
properties of these spaces. 


For v > — 1 let A v : C —> C and B v : C —> C be real entire functions defined by 

00 (-1 ) n (\ Z ) 2n 


A A z ) = 


n =0 


i\{v + l)(v + 2)... (i ' + n) 


and 


B A Z ) = 




\ 2n+l 


n—0 


n\(v + 1 )(y + 2)... (v + n + 1)' 


These functions are related to the classical Bessel functions by the identities 


A v {z) = + 1) {\z) " J„{z), 

B„(z) = T(v + 1) (iz) v J v+ \(z). 


( 1 . 11 ) 


( 1 . 12 ) 


(1.13) 

(1.14) 


It follows that both A v and B v have only real, simple zeros and have no common zeros (note that in the 
simplest case v = —1/2 we have A_ x / 2 {z) = cosz and B_ 1 / 2 (z) = sinz). Moreover, they satisfy the system 
of differential equations 


A' u (z) =-B„(z) and B/(z) = A u (z) - (2v + l)z 1 B v (z). 


(1.15) 


The entire function 

E v (z) := A u (z) - iB v (z) 

is a Hermite-Biehler function that satisfies properties (PI) - (P4) listed above, with v{E v ) = t^E^) = 1 (we 
shall verify these assertions later), and we denote by K v the reproducing kernel of the de Branges space 
U(E V ). 

For S > 0 and N £ N we define to be the set of all entire functions F : C N —> C of exponential type 
at most S. Given a real-valued function T : —y K. we write 


Ef-(J-) := {C £ Ef; £(x) < .F(x), Vx £ 1^}, 
E*+(F) ■.= {M £ Ef; M(x) > F(x),Vx £ R N }. 


Theorem 2. Let v > —1, A > 0 and F\(x) = e 7rA l x l 2 . Let 

t/^-(<5,A) = inf|^| e -^W 2 -£(x)} |x| 2 ^ 2 - w dx; C £ Ef"(J\)| , 

U^ + (6, A) = inf | j |a4(x) _ e -^l x l 2 J | x | 2 ^+ 2 -iv dx; m £ Ef^(J’ a) 

The following properties hold: 

(i) If k > 0 then Uf r± (S, A) = k 2i ' +2 Uff ± (n5, n 2 X). 

(ii) For any S > 0 we have 

Uu ± (S, A) = 5 ujn-i Ul ± (8, A), 

where ojn -i = 2tt n ^ 2 T(N/2)~ 1 is the surface area of the unit sphere in . 
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(iii) When N = 1 and 5 = 2 we have 




Ui~{ 2,A) = 

r(^ + i) y- 

e —7 tA|^| 2 

(1.16) 

(7TAU+ 1 ^ 

c*,AV (£,£)’ 

Ul + { 2,A) = 

v e- A lf 2 

i> + i) 

(1.17) 

B „«)= 0 CvKv ^'0 

(7rA) l/+1 ’ 

w/iere c v = -k 2 _2i,_1 T(v + l) -2 . 





(iv) There exists a pair of real entire functions of N complex variables, z i— > C V>N (S, A,z)£Ef-(J*) and 
z i—^ A, z) £ Ef+pV), which are radial when restricted to and verify 

U?~(6, A) = [ {e"^l*l a -Z: i/>J v( ( y ) A ) x)}| X | 2 ^ 2 - Jv dx 

and 

U* + (6,\) = [ {a^^jv(^A,x)- e _7rA|x|2 | |x| 2y+2_Jv dx. 

J R w L > 


In order to extend the construction of Theorem[5]to a class of radial functions, it is important to understand 
the asymptotics of (11.1611 and (11.171) as A — > 0 and as A — > +oo. From (11.71) . (11.131) . (11.141) and (11.151) we 
have, for £ >0, 

= 2- 1 c” 1 r 2!/_1 {^(0 2 + C^+i(0 2 ~(2p + 1) MO Mid)}- (1-18) 

Using the well-known asymptotic for Bessel functions 

MO = cos (£ - ¥ - f) + 0,(r 3/2 ) (1-19) 

as £ —> oo, we have 

hm {^(^) 2 + £J„+i(0 2 -(2^ + 1) J v (OMi (£)} = f (1-20) 

It then follows that 

K v (t,Q~n- 1 c; 1 r 2v - 1 ( 1 - 21 ) 


as £ —>■ oo. Also, notice by (11.131) and (11.141) that the positive zeros of A v and B u are just the positive zeros 
of J u and Jy+i, and it is well-known that, for large m, the Bessel function J v has exactly m zeros in the 
interval (0,m7r + ^ + f ] (see [Ml Section 15.4]). Denoting by j„,n the n-th positive zero of J v (ordered by 
size), it follows from (11.161) and (11.211) that 

i 00 

|#"(2,A)| «„ ^ + E nMle ' TA, ''“ 

n— 1 

i 


as A —> oo (we use Vinogradov’s notation / <C j, or / = O(g), to mean that |/| is less than or equal to a 
constant times \g\, and indicate the parameters of dependence of the constant in the subscript). Analogously, 
from (11.171) and (11.211) we find 


M + ( 2, A) I «„ 


1 

c„K v { 0,0) 


+ n2 " +1 e~" Xj ’ +1 ’ n 

n= 1 


i 


+ 


l 

ATT 1 
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as A —> oo. To complete the asymptotic description we shall prove in the beginning of Section [5] that 

\U^( 2,A)|<„, fc A fc (1.22) 

as A —> 0, for every k G N. It is precisely this fact that makes the Gaussian subordination framework more 
robust (when compared to the analogous framework generated by exponential x i —> e _A l x in [TT] . that verifies 
(11.221) only for k = 1). A somewhat surprising feature is the fact that the solution for the exponential is an 
important step towards the solution for the Gaussian, as we shall see in Section [2j 

1.5. Gaussian subordination. We now proceed with the extension of Theorem [2] to a class of radial 
functions. The way we define our class of functions below might look slightly odd at first, but it will be clear 
later on that this is indeed the most general way to proceed. This is essentially motivated by the fact that, 
in order to make the best use of the strong decay at the origin given by (11 .22|1 . we are led to work on the 
Fourier transform side. In general, this framework covers functions that are not integrable. This requires us 
to work with Fourier transforms in the distributional sense. 

We consider two classes of nonnegative Borel measures g on (0,oo) for a given parameter v > — 1. For 
the minorant problem we consider the class of measures /i satisfying the condition 

/•OO 

J o 1 + Xv+k+1 < oo (1.23) 

for some fceN, whereas for the majorant problem we consider the class of measures p satisfying the more 
restrictive condition 

roc \k 

l 1 + A fc d^(A) < oo (1.24) 

for some k € N. Let g^ : R —> R U {oo} be a function satisfying the following properties: 

(Ql) The function g^ is even. 

(Q2) For the majorant problem, the function g M is finite-valued and continuous on R. For the minorant 
problem, the function g^ is finite-valued and continuous on R \ {0} and g fl ( 0) = limsup x _ >0 g^(x). 

(Q3) For each N £ N, let G/j.,n '■ —► R be the radial extension of g M to R^, i.e. 

G/i,n(x) = fl M (|x|)- (1-25) 

If N v = \2v + 2] (the smallest integer greater than or equal to 2v + 2), then G^,n u is a tempered 
distribution on M' Vl/ . 

(Q4) For any Schwartz function <p : R iV “ —y C that is zero in an open neighborhood of the origin we have 
[ Gii,n u (x) ^(x)dx = [ (f \~ N »/ 2 e — 7r l t l 2 / A d/x(A)'] ^(t) dt, (1-26) 

JJb. n » \J o J 

which, informally speaking, means that we require the Fourier transform of Gh,n„ to equal the 
expression in parentheses in (11.261) for t ^ 0. 

For K>0we define a new measure /i K on the Borel subsets X C (0, oo) by putting 

g K (X) := n( K X), (1.27) 

where 

kX = {kx\ x € X}. 

Note that n K satisfies (11.231) or (11.241) whenever n does. We are now in position to state our next result. 
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Theorem 3. Let p be a nonnegative Borel measure on (0,oo) satisfying (11.231) for the minorant problem or 
(11.241) for the majorant problem. Let gy :1->KU {oo} be a function satisfying properties (Ql) - (Q4) above 
and let G^^n '■ R N —> R be defined by (11.251) . For N £ N write 

U?~(5,p) = inf {^ {^,jv(x) - £(x)} \x\ 2v+2 ~ N dx; £ £ K~(G^n) j , 

U? + (8,p) = inf [J n |A4(x) - at(x)| Ixl 2 ^ 2 ^^ dx; M £ Ef + (^ jAr )| . 

The following properties hold: 

(i) If 0 < k then U^ ± {8,p) = n 2v+2 U^^(k 8, p K -f)- 

(ii) For any S > 0 we have 

p OO p OO 

U^(8,p) = U^(S, A)d/i(A) = |w JV _ 1 / Ul ± {8, A) dp(X). (1.28) 

Jo Jo 

where u>n-i = 2n N ^ 2 T(N/2)~ 1 is the surface area of the unit sphere in M. N . 

(iii) There exists a pair of real entire functions of N complex variables, z i— > £v,n(8,h, z) £ E^ (G/j.,n) 
and z i-4 M u ,n{ 8, p, z) £ E {G^,n), which are radial when restricted to R^ and verify 

U*~{5,p) = [ {^.jv(x) — £ v , N ( 8, p, x)} |x| 2,/+2_jv dx 

J 

and 

U? + {8,p) = [ p, x) — G^nOO} \x\ 2l, ^ 2 ~ N dx. 

Remark: When restricted to the class of entire functions of N complex variables that are radial on R^, the 
extremal functions described in Theorem [2] (iv) and Theorem [3] (iii) are unique. This follows from the recent 
work of F. Gongalves m Section 4.2] on interpolation formulas with derivatives in de Branges spaces. On 
the other hand, uniqueness of the extremal functions described in Theorem[l]is only known under additional 
conditions on E (see j2T] Theorem 1]). 

The idea of solving the extremal problem for a base case with a free parameter A > 0, and then integrating 
the parameter against a suitable measure p is reminiscent of the work of Graham and Vaaler I23j . As the 
final result (11.281) suggests, a necessary condition to be imposed on the measure p is the /i-integrability of 
Ul ± {8, A) considered as a function of A. In order to generate the largest possible class, one wants to make 
this necessary condition also sufficient, and that is precisely the reason why we define our objects via the 
Fourier transform (since the Fourier transform of the Gaussian in (11.261) has an exponential decay to zero as 
A -4 0 for each fixed t ^ 0). The passage to dimension N v in (11.261) guarantees the proper decay as A —► oo. 
In Section [6] we describe a variety of examples that can be obtained via our Gaussian subordination method. 
In particular, this new framework strictly contains the class m of our precursor m- 

The one-dimensional Gaussian subordination framework b m provided the solution of the Beurling- 
Selberg extremal problems for previously inaccessible functions such as f a ,b(x) = log ((a; 2 + b 2 )/(x 2 + a 2 )), 
where 0 < a < b; g(x) = arctan(l/x) — x/(l +x 2 ) and h(x) = 1 — a;arctan(l/a;). These were used to improve 
the existing bounds for the modulus and the argument of the Riemann zeta-function on the critical strip 
under RH (also for L-functions under GRH) [4] j6 f [7]. We expect that the extension of this framework to a 
multidimensional setting presented here might reveal other interesting applications. 














2. Entire functions with prescribed nodes 


2.1. Laguerre-Polya functions. We recall that a Laguerre-Polya function F is a real entire function with 
Hadamard factorization given by 


F{z) ^Em z r e -^ 


oo 




( 2 . 1 ) 


3 =1 


where r £ Z + , a, &, £ K, with a > 0, ^ 0 and Y^jLi < 00 (the product may be finite or empty). We 

say that F has degree Af = Af(F), with 0 < Af < oo, if a = 0 in (12.11) and F has exactly AT zeros counted 
with multiplicity. Otherwise we set Af(F) = oo. 


The important property of Laguerre-Polya functions for our purposes is that we can write their inverses 
as Laplace transforms on vertical strips. In fact, if ( 71 , 72 ) CM is an open interval not containing any zeros 
of F (we allow ti,T 2 £ {±oo}) and c £ (ti,T 2 ), we have 


F(z) 


f 


g c (t) e zt d t 


for 7”i < Re (z) < r 2 , where 


1 f 

9c{t) = ^l 


c+ioo st 


F{ s ) 


d s. 


( 2 . 2 ) 


(2.3) 


The integral (12.31) is absolutely convergent if Af(F) > 2, and is understood as a Cauchy principal value if 
Af(F) = 1. If Af(F) = 0, (12.21) holds with g c being an appropriate Dirac delta distribution (though we will 
not be interested in this case here). A simple application of the residue theorem gives us that g c = gd if 
c, d £ ( 71 , 72 ). We briefly review below the main qualitative properties of the frequency functions g c needed 
for this work. We refer the reader to [26] Chapters II to V] for a detailed account of this theory. 


Lemma 4. Let F be a Laguerre-Polya function with at least one zero and let g c be defined by (12.31) . where 
c £ K and F{c) ^ 0. The following propositions hold: 

(i) The function g c is real-valued and of one sign, and its sign is equal to the sign of F{c). 

(ii) If F has degree Af > 2, then g c £ C'^“ 2 (IR.). 

(iii) If ( 71 , 72 ) C R is the largest open interval containing no zeros of F such that c £ ( 71 , 72 ), then for 
any t £ ( 71 , 72 ) we have the estimate 

|< 7 c(t)| -Cr e Tt Vf £ M. (2.4) 

Proof. Part (i) follows from [2S] Chapter IV, Theorem 5.1]. For part (ii), observe that we can differentiate 
under the integral sign in (12.31) since 

l/\F(c+iy)\=0(\y\~ k ) (2.5) 

as \y\ —> 00 if Af(F) > k (for details see [2|j] Chapter II, Theorem 6.3 and Chapter III, Theorem 5.3]). For 
part (iii), observe that g c = g T and estimate the integral in (12.31) trivially using (12.51) (in the case Af{F) > 2). 
If Af(F) = 1, we can see it directly as in (12.91) below. □ 


We use the notation g c _ for the function g c _ e in (12.31) . with £ so small that there is no zero of F in the 
interval (c — e, c), and analogously for g c+ . In this paper we are especially interested in the case when the 
nonzero roots {A ; }°Li of F verify 

OO 

> 0 for all j and E f- 1 < (X). (2.6) 

3 =1 
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The next lemma presents an essential property of the Laguerre-Polya functions satisfying m- 


Lemma 5. Let k G {0,1}. Let F be a Laguerre-Polya function of the form 

F(z)=z k f[(l--f), (2.7) 

i=i V J 

where the points satisfy (12.61) . Then the function go- defined in (12.311 satisfies go-(t ) = 0 for all 

t > 0. In particular, m becomes, for Re (z) < 0, 

<2,8) 


Proof. 


Let us consider first the case k = 0. Assume that - and define 


*■-(*>=n (! - j- 


For a > 0, we define the function y a _ by 


7a- (t) = 


ae at , if t < 0; 


0, 


if t > 0. 


We note that 


, — = I 7 a-{t)e zt dt 

1 a J - oo 

in the half-plane Re ( 2 ) < a. It is then straightforward to check that 

1 f°° 

rrr = / ( 7 ^- * • • • * 7 £„-)(*) e zt dt 

n\Z) J— OO 


F n {z) 


(2.9) 


in the half-plane Re (z) < £ 1 . Note that all functions t >->■ 7 a _(t), for a > 0, belong to L 1 (R) D L°°(R), 
hence the convolution t > ( 7 ^_ * ... * 7 £„ — ){t) is a well-defined continuous and integrable function for 
n > 2. Since 1 /F n converges uniformly to 1/F on the line Re (z) = 0 and > 0, it follows from (12.31) that 
( 7 £j_ * ... * 7 ^ n _) converges pointwise to go as n —> 00 . In particular, go(t) = 0 for t > 0. Since go = go-, 
the result is shown if k = 0. If k = 1, the n-fold convolution is convolved with the function that equals — 1 
for t < 0 and 0 for t > 0. We use in this case that 1 /F n converges uniformly to 1/F on the line Re (z) = —£. 
The calculations are analogous. □ 


2.2. Interpolating the exponential. Let F be a Laguerre-Polya function given by (12.71) where the points 
satisfy (12.61) . Let g = go- as in (12.31) and define, for A > 0, the functions 

A(F, A, z) = e~ Xz - F(z) f g(w - A) e~ zw d w, for zeC, 

0 J ° ( 2 - 10 ) 

A 1 (F,X,z) = F(z) g(w - A) e~ zw dw, forRe(^)<0. 

J —OO 


Since the smallest zero of F is nonnegative, Lemma @] (iii) implies that z 1 —> A(F, A, z) is entire and z 1 —>• 
Ai(F, A, z) is analytic in the half-plane Re (z) < 0. Also, it follows directly from (12.81) that 


A{F, A, z) = A\(F, A, z) 


( 2 . 11 ) 


in Re (z) < 0. 
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We now define, for z € C and A > 0, the functions 


H(z) = 


I(F,\,z) = 


1, if Re (z) > —1/4; 
0, if Re (z) < -1/4; 
r \ 


-f g(w- X)e~ zw dw, if Re(z) > -1/4; 
Jo 

[ g(w- A) e~ zw d w, if Re (z) < -1 /4. 

J — OO 


The value —1/4 is not essential in the definitions above; any r < 0 could be used in its place. Equations 
(12.101) and (12.111) can be rewritten as 

A(F, A, z) - H(z)e~ Xz = F(z)I(F, A, z) (2.12) 

for all z £ C. We establish in the next two lemmas some straightforward but crucial properties of these 
functions. 


Lemma 6. The functions z K>• I(F, A, z) and z K>• H(z)e Xz are bounded in the complex plane (by a constant 
depending on A). 


Proof. We consider the function I. For Re (z) > —1/4 we have 


\I{F, A, z)\ = 


[ g(w - A) e~ zw dw 

Jo 


Oa(1). 


It follows from Lemma [S] that g = go- is identically zero for t > 0, and it follows from Lemma [I] (iii) (with 
r = —1/8, say) that \I(F, A, z)\ < C|Re(^) + 1/8| —1 for Re (z) < —1/4. The claim for z i-> H(z)e~ Xz is 
evident. □ 


Lemma 7. Let £ = ±1 be the sign of F in the interval (—oo,0). For all x £l we have 

eF(x){e~ Xx -.4(F,A,x)} > 0. (2.13) 

Proof. From Lemma H (i), the function g is of one sign, equal to the sign of F in (—oo,0) (in fact, a 
considerably stronger statement regarding the derivatives of g is true, cf. Chapter IV, Theorems 5.1 
and 5.3], but we do not require it here). The inequality now follows directly from the definition of A in 

dzm □ 


Note that, with our normalization (12.71) . we have e = 1 if -F(O) ^ 0 and e = —1 if F( 0) = 0. 


2.3. Interpolating the Gaussian. Let A : C —> C be an even Laguerre-Polya function of exponential type 
t(A). Assume that A(0) = 1 and label the positive zeros of A by 

0 < ai < a 2 < .. . 

with repetition according to multiplicity. Note that for all z € C the Hadamard product representation 










is valid (the absence of exponential factors in the product is crucial for our purposes). The product converges 
uniformly and absolutely in compact sets of C. Define F A : C —> C by 


z 


(2.14) 


f A {z) = n 1 - 

1=1 V “■?/ 

and observe that F A is also a Laguerre-Polya function, and that it verifies 

A(z) = F a (z 2 ) 

for all z € <C. 

In a similar way, let B : C —> C be an even Laguerre-Polya function of exponential type t(B), with a 
double zero at the origin, and positive zeros listed as 


0 < 6i < fc 2 < • • • 


If we assume that B"( 0) = 2, we have 


b(z)= z 2 n 1 - 

i=i V 


ST ■ 


and we may define Fb : C —1 C analogously to (12.141 ) in order to have B(z ) = Fb ( z 2 ). Observe that F A and 
Fb verify (12.61) and (12.71) : F A with k = 0 and Fb with k = 1. 


Lemma 8. Define functions L and M by 

L(A, A, z) = A(F Ai tt\,z 2 ); 

M(B,A,z) = A(F b ,ttA,z 2 ). 

The following propositions hold: 

(i) The function z L(A , A, z ) is an entire function of exponential type at most t(A), and for every 

x € R. the inequality 

A(x){e~ nXx2 - L(A, A,x)}>0 (2.15) 

holds. Moreover, for all f with Al(£) = 0 we have 

L(A,X,f) = e~^ 2 . (2.16) 

(ii) The function z <—y M(B,A,z) is an entire function of exponential type at most t(B), and for every 
x £ M the inequality 


B(x){m(B,\,x) -e-^ 2 } 

holds. Moreover, for all f with B(£) =0 we have 

M(B, A, £) = e~^ 2 . 


> 0 


(2.17) 

(2.18) 


Proof. Since 2 A{F a , 7tA, z) is entire, it follows that z L(A, A, z) is entire. From (12.121) we have 

L(A, A, z) = H (z 2 ) e ~^ z2 + F a (z 2 ) I(F a , ttA, z 2 ), 

and Lemma [6] implies that there exist constants C,D > 0 such that 

\L(A,X,z)\<C + D\F a (z 2 )\=C + D\A(z)\ 
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for all z?C. It follows that z L(A, A, z ) has exponential type at most t(A). The inequality (12.1511 follows 
from Lemma [TJ while identity (12.161) follows from the definition of A in (12.101) . This completes the proof of 
(i). The proof of (ii) is analogous. □ 


3. Proof of Theorem [0 

The general strategy we follow here is similar to [11, Section 3]. Let £ be a Hermite-Biehler function 
satisfying properties (PI) - (P4), and let A and B be the companion functions defined by (11.61) . Property 
(PI) implies that A and B are Laguerre-Polya functions of exponential type at most t(E). Property (P3) 
implies that A is even and B is odd, while property (P2) implies that A(O) ^ 0 and B has a simple zero at 
the origin. Further details are given in m Lemmas 12 and 13]. 

Lemma 9. Let A > 0. Let z L(A 2 , A, z) and z >->• M(B 2 , X,z) be defined as in Lemma[8l Then 

x !->■ L(A 2 , A, x ) — e~^ Xx 

and 

x^M(B 2 , \,x)-e~ nXx2 

belong to L 1 (M, \E(x)\~ 2 dx). 


Proof. We note from (12.101) that 


since g is bounded. This implies the first claim since \A/E\ is bounded by 1 on the real line. The second 
claim is shown analogously. □ 


L(A 2 , A, x) — e 


— TVXX 2 


A 2 {x) 


g(w — 7tA) e 


"die 


Lemma 10. Let E be a Hermite-Biehler function of bounded type in U with exponential type t(E). 

(i) If M : C —> C is a real entire function of exponential type at most 2 t(E), which is nonnegative on R 
and belongs to Ir(R, \E(x)\~ 2 dx), then there exists W £ TL(E) such that 


M{z) = W{z)W*{z) 


(3.1) 


for all z £ C. 


(ii) In particular, if A > 0 and 



|.E(ir)| 2 dx < oo, 


there exist U, V £ H{E) so that 


M(B 2 ,X,z) = U(z)U*(z), 

L(A 2 , A, z) = U(z)U*(z) - V(z)V*(z) 


for all z £ C. 


(3.2) 


(3.3) 

(3.4) 


Proof. Part (i) is m Lemma 14]. To prove part (ii), note by Lemma [8] that z i-)- L(A 2 ,A, z) and z >->• 
M(B 2 , A, z) have exponential type at most 2 t(E). From Lemma|9]and (13.21) . it follows that x i —> M(B 2 , A, x) 
is an element of L 1 )®., |£’(a:)| _2 da;) and this gives the representation (13.31) as a particular case of (El- 
Representation (13.41) follows similarly, considering the real entire function z i-» M(B 2 , X, z) — L(A 2 1 X,z), 
which is nonnegative on R. □ 
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Proof of Theorem QJ Let M : C — > C be an entire function of exponential type at most 2 t(E) such that 
M(x) > e-™ Xx2 for all x £ R and 

/ OO 

M(x) \E(x)\~ 2 da: < oo. 

- OO 

Lemma flOl implies that M = WW* for some W £ TL(E). Since B ^ TL{E ) by hypothesis (P4), we are in 
position to use [3 : , Theorem 22] to conclude that {z K{^,z)\ B(£) = 0} is a complete orthogonal set in 
TL(E). By Plancherel’s identity we have 


f 


M(x) |i?(a;)| 2 dx = 


pOO 

W(x) 

} —oo 

E{x) 


dx= 


B (5)=0 


l^(Q | 2 


> E 


B(0 =0 


e -n\i 

Wx' 


(3.5) 


which proves (11.101) . From Lemma|S]we know that z i—>■ M(B 2 , A, z) is an entire function of exponential type 
at most 2 t(E) that satisfies 


M(B 2 , \,x) > e~* Xx2 


for all x £ R. Moreover, it follows from (12.181) and (13.31) that we have equality in (l3~5l) for M(z) = M(B 2 , A, z). 


Now let L : C —> C be an entire function of exponential type at most 2 t(E) such that L(x) < e nXx for 
all x £ R and 


/ OO 

L(x) |-E(x)|~ 2 da; < oo. 

-OO 


We use the fact that there exists a majorant z M(B 2 ,X,z) in L 1 (K, |£?(a:)| -2 da;). This allows us to 
apply Lemma [10] to the function z i-» M(B 2 ,X,z) — L(z ) (which is nonnegative on the real axis, belongs 
to L 1 (R, |£'(x)|" 2 da;) and has exponential type at most 2 t(E)). Together with (13.31) . this shows that L = 
UU* — TT* with U,T £ TL{E). Since A TL{E) by hypothesis (P4), a new application of [3J Theorem 22] 
shows that {z K > K(£,z); A(£) = 0} is also a complete orthogonal set in TL(E). Plancherel’s identity then 
gives us 


f 

J —( 


L{x)\E{x)\ 2 dx = 


f 

J —( 


\U(x)\ 2 ~\T{x)\ 2 
\E(x )\ 2 


d x = yy 
> 1 ( 0=0 


TOi 2 -roi s 

mx 


s E 


A(€)=o 


KXX 


(3.6) 


which proves m- From Lemma [5] we know that z K > L(A 2 , A, z) is an entire function of exponential type 
at most 2 t(E) that satisfies 

L{A 2 , X,x) < e~* Xx2 

for all x £ R, and it follows from (12.161) and (13.41) that we have equality in (13.61) for L(z) = L(A 2 , A, z). □ 


4. Proof of Theorem [2] 

4.1. Preliminaries. 

4.1.1. Homogeneous spaces. We start by recalling the basic properties of the spaces TL(E^) : which are ex¬ 
amples of homogeneous de Branges spaces j3] Section 50]. Further details are given in [27l Section 5] and 
[TT] Section 4]. 

Recall that v > — 1 is a given parameter and 

E v {z) := A v {z) - iB v (z), 

with A v and B u given by (11.111) and (11.121) . respectively. From [3J Problems 227 and 228] we have that E v 
is in fact a Hermite-Bichler function with no zeros on the real line. Moreover, E v has bounded type in IA 
with v(E v ) = t(E v ) = 1. By definition, A v is even while B u is odd. This plainly implies that z K > E^iiz) is 
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real entire. This accounts for properties (PI) - (P3) for E v . The next lemma collects the relevant properties 
of the spaces H{E U ) for our purposes. 


Lemma 11. Let v > — 1. The following properties hold: 

(i) There exist positive constants a„,b u such that 

a v \x\ 2v+1 < | E v {x)\~ 2 < b v \x\ 2 " +1 (4.1) 

for all x € R with |a:| > 1. 

(ii) For F £ "H(Ify) we have the identity 

/ oo roo 

\F(x)\ 2 \E„(x)\~ 2 dx = c v / \F(x)\ 2 |x| 2ly+1 dx , (4.2) 

-oo J — OO 

with c v = 7T 2 _2 " _1 r(i/ + lfy 2 . 

(iii) An entire function F belongs to R(E V ) if and only if F has exponential type at most 1 and 

/ OO 

|F(a;)| 2 |ir| 2l ' +1 da; < oo. (4.3) 

-OO 

Proof. We refer the reader to E3 Section 5]. □ 


In particular, the identity (14.211 is fundamental for our purposes, and provides a good example of the 
power (and a bit of mystery) of the de Branges space machinery. From (11.131) , (11.141) and (11.191) we see that 
A v and B„ do not verify (|4.3I) and thus A V ,B„ £ R{E U ). This establishes the remaining property (P4). 


4.1.2. Symmetrization techniques. Before we move on to the proof of Theorem [2] we also to need to recall 
the basic symmetrization mechanisms to deal with radial functions in several variables. This is contained in 
m Section 6] and we again replicate here a brief collection of the main facts for the reader’s convenience. 

We start with our extension operator. If F : C —> C is an even entire function with power series 
representation 

OO 

F(z) = J2 ^kz 2k , 
k=0 

we define the extension : C N —> C by 

oo 

ip N (F)(z) = c fc(~i + • ■ • + z 2 N ) k . 
fc=o 

We also consider a suitable symmetrization operator. For an entire function J- : C N —> C, with IV > 1, we 
define its radial symmetrization T : C N —> C by 

E(z) = j F(Rz)da(R), 

JsO(N) 

where SO(N) denotes the compact topological group of real orthogonal N x N matrices with determinant 1, 
with associated Haar measure a normalized so that a(SO(N)) = 1. If IV = 1 we set E(z) = i{J r (z)+J r (—.j)}. 

Lemma 12. The following propositions hold: 

(i) Let F : C —> C be an even entire function. Then F has exponential type if and only if if n{F) has 
exponential type, and t(F) = r(fyv(F)). 
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(ii) Let T : C N C be an entire function. Then T : C N —> C is an entire function with power series 
expansion of the form 

OO 

=J2 c k(zi+... +z%) k . (4.4) 

k =0 

Moreover, if J- has exponential type then J- has exponential type and t(.F) < t(JF). 

Proof. These are E3 Lemmas 18 and 19]. □ 

4.2. Proof of Theorem [2]. Recall that we write .F A (x) = e _7rA l x l 2 . 

4.2.1. Proof of part (i). For k > 0 observe that C 6 E^ _ (J r A) if and only if x t->- £(kx) 6 E^ _ (J' k 2 A ). An 
analogous property holds for the majorants. This is enough to conclude part (i). 


x\ 2v+1 dx < oo. 


4.2.2. Proof of part (ii). Let LgEj (,F A ) be such that 

/ °° . . 

{e-^ 2 -L(x)\ 

-OO ^ J 

By replacing L by L , we may assume that L is even, and hence has a power series expansion 

OO 

L{z) = y ^C k z 2k . 

k—0 

By Lemma [T^l fil we have that i/jn(L) G E^ _ (.F a ), and a change to radial variables gives us 

J {e ~^ 2 -^(L)(x)} |x|^+ 2 -^dx = I™ {e ~” Xx2 - L(x)} {x ] 2 ^ 1 dx. 


Hence 


U?~{8, A) < \u N -xUl~{8, A). 


On the other hand, let £ G E^ (T\) be such that 

[ | e —irA|x| 2 _ £( x )| | x |2i/+2-iV dx < OO. 

By Lemma fl2l (ii) we have that C G E^ _ (.Fa) and that it has a power series expansion 


£( z ) = y] c fc( x i +--- + 4) fc 

Define the entire function L : C —> C by 


k—0 


L(z) = C{z, 0,0,...,0) = ^ 


C fc z 2fe . 


fc =0 


Another application of Lemma [T2] (i) gives us that L G Ej (fP\) and we arrive at 


cov-i r \e~* Xx2 ~L{x)}\x\ 2, ' +1 dx= [ |e _7rA l x l 2 

J-oo L J Jr n L 



— £(x) j Ixl^+^dx 

— £(x) j Ixl^-^dx. 


Hence 

iu N . 1 Ul-(5,X)<U^-(S,X), 

and this concludes the proof for the minorant case. The majorant case is treated analogously. 
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4.2.3. Proof of part (iii). We are now interested in computing U„(2,A). Let us start with the majorant 
case. First observe that 

— 2 |d 2 - 1 dx= I >± 1 ) 


/ 

J —( 


(7tA) 1/+1 ' 

Let M : C —> C be an entire function of exponential type at most 2 such that M[x) > e _7rAx for all x £ 
and 

/ OO 

M(x) \x\ 2l/+1 dx < oo. 

-OO 


From (14.11) we obtain 


/ OO 

M (x) \E v {x)\~ 2 dx < oo. 

-OO 


It follows from Lemma [TU] that M = WW* for some W £ 'H(E V ). From Theorem [T] and the key identity 
(IQ) we have 

0 —77A{ 2 


b„(€)=o 


/■ oo 

KAU) ~ L M{x) lK{x)r2 dx 

\W{x)\ 2 \E„{x)\- 2 dx 

/ OO 

\W(x)\ 

-oo 

M(x) 


(4.5) 


2 |x| 2l/+1 dx 


' —OO 
poo 


|x| 2 " +1 dx. 


Moreover, we have seen that the entire function z M(U^, A, z) of exponential type at most 2 realizes the 
equality in (14.51) . Therefore 

e"^ 2 I> + 1) 


Z7„ 1+ (2,A)= ]T 


c v K v {^,i) (ttA)^ 1 


0 

The minorant case follows analogously, using the fact that there exists a majorant z K > A, z). This 

allows us to write a general minorant L as UU* — TT* : as done in the proof of Theorem [T] 


4.2.4. Proof of part (iv). For N > 1 and 5 = 2 we define £,,^(2, A, z) := iI>n(L(AI, A, •)) (z) and A, z) 

A, •)) (z). These functions have exponential type at most 2 and satisfy 

£„iv(2,A,x) < e _ 7 rA l x l 2 < Mv , n { 2,A,x) 

for all x £ K ;V . Moreover, 

J n |Al^]v(2, A,x) -e _,rA|x|2 | |x| 2l ' +2 ~ Ar dx = \u N -i J |m(B 2 ,A,x) - e -7 ^ 2 } |x| 2l/+1 dx 

= |wjv-if/ ! ) + (2,A) 

= U?+(2,\). 

The computation for the minorant is analogous. The existence of extremal functions for general S > 0 follows 
by the change of variables given in part (i). This completes the proof. 
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5. Proof of Theorem 0- Gaussian subordination 


5.1. Preliminaries. We start by showing that the asymptotic estimate (11.221) holds. We have established 
the integral representation 

/ OO j pTT \ | 

l Bl{x) J g{w — 7 rA) e~ x2w dw > |x| 2 l,+ 1 dx, 

where g = go- for the Laguerre-Polya function F B 2 as in (12.31) and (|2.14D . We note that B 2 /\E V \ 2 is bounded 
by 1 on the real line, and therefore, using (ED, we see that B 2 (x) \x\ 2u+1 is also bounded on the real line. 
Moreover, g(t) = 0 for t > 0 and g £ C' 00 (R). This implies, in particular, that g^ k \ 0) = 0 for all k £ N, 
which leads to | g(w — 7rA)| < Cfc(7rA — w) k for A < 1, say, and all 0 < w < ttX. An application of Fubini’s 
theorem gives (11.221) for the majorant. The minorant case follows analogously using A„. 

If the nonnegative Borel measure g on (0,oo) satisfies (11.231) we define 

roc 

= J |e _7rA|x|2 - C v ,n{6, A,x)j dg{X), (5.1) 

and if g satisfies (11.241) we define 

roo 

= J {^.N(^A,x)-e-^ |x| -}d / r(A), (5.2) 

where z 1 —> £^^(<5, A, z) and z 1 —> A, z) are the extremal functions of exponential type at most <5 de¬ 

fined in Theorem[2](iv). Integration in x with respect to |x| 2l/+2_Jv dx and an application of Tonelli’s theorem 
implies that the nonnegative functions x D± N (5, g, x) are radial and belong to L 1 (R JV , |x| 2l/+2 ^ JV dx), 
hence the integrals m and m converge almost everywhere. 

Lemma 13. The following properties hold: 

(i) If /i satisfies (11.231) . then the nonnegative function x 1 —» D~ N (S, g,x) is continuous for x ^ 0. 
Moreover, D~(5,g,x) = 0 whenever ^4^(<5|x|/2) = 0. 

(ii) If g satisfies (11.241) . then the nonnegative function x 1 —> D+ N (5, g, x) is continuous. Moreover, 
D+ (5, g,x) = 0 whenever B u {5 |x|/ 2 ) = 0 . 

Proof. We restrict ourselves to the case <5 = 2. The general case follows by a standard dilation argument. 

Part (i). Recall that 0) ^ 0. Let g = go- for the Laguerre-Polya function F 42 as in (12.31) and (12.141) . 
and note by Lemma []] that g > 0 on R. We have 

/* 7T A 

e -7i-Ax — L(A 2 , A, x) = A 2 (x) / g(w — nX) e~ x w dw. 

J 0 

Since g > 0 on R it follows that Hfi defined by 

p7rX 

Hf(x)= / g(w — nX) e~ x w dw 

J 0 

is nondecreasing on (— 00 ,0] for every A > 0. For x < 0 we define 

pOO 

df(2,g,x):=Al(x) Hfi(x)dg(X), 

Jo 

and note that this function is continuous on (— 00 , 0 ) by dominated convergence (we may have 
limsup x ^ 0 - d~(2, g, x) = 00 ). It follows that df(2,g,f) = 0 for every £ < 0 with A„(f) = 0. Since 
^v,N ( 2 ,/i, x) = d v ( 2 ,/i, —|x|), our claim follows. 
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Part (ii). Similarly, consider g = (?o_ for the Laguerre-Polya function F B 2 as in (12.31) and (12.141) . Now we 
have g < 0 on R and 

PTC X 

M(Al 2 , A, x) — e~ nXx =—Bl(x) / g(u> — 7 tA) e~ x w dw. 

Jo 

Define 

PTC X 

H^ (x) = — / g(w; — 7 rA) e~ x w dw 

Jo 

and note that this is a nonnegative function which is nondecreasing on (—oo,0] for every A > 0. For x < 0 
define 

p OO 

dt(2,n,x) := Bl(x) / i?+(x)d^(A), 

Jo 

and again, we note that this function is continuous on (—oo, 0) by dominated convergence. Thus (2, g, £) = 
0 for every £ < 0 with B„(t;) = 0 . 


To analyze the situation as x — > 0 , we use integration by parts to get 

ff(- ttA) 1 




fK A 

■ / g'(w — 7 tA) e” 

Jo 


’ dw, 


and thus 


d+(2,g,x) = B ^f > J g'{w-irX)e x w dw j dg(X). 


Observe now that 


and 


\g(-Tr\)\dg(X) < oo 


lo 


(5.3) 

(5.4) 

(5.5) 


P OO PTC A 

/ / |(/(u; — 7rA)| dw dg(X) < oo. 

Jo Jo 

In fact, (15.41) follows from the fact that g is bounded and g^ (0) = 0 for any k > 0. To arrive at (15.51) . we use 
the additional fact that |c/'(t)| -C e Tt for some r > 0 (Lemma U (hi)), which holds since g' is the frequency 
function associated to the Laguerre-Polya function F B 2 ( z ) /z on a vertical strip containing the origin. We 
remark that F b i (z) /z does not have a zero at the origin anymore. We may then apply dominated convergence 
in (15.31) to conclude that 


lim d+(2,g,x) = { lira — ^ 1 [ (—g{— 7 rA) — f g'(w — ttX) dw | dg(X) 

x^-o- x z ) J Q y J 0 J 


= 0 . 


The result now follows since D+ N (2, g, x ) = d+(2, M ,-|x|). 


□ 


The following lemma is crucial to prove the optimality part of Theorem [3] 

Lemma 14. Let v > — 1 and G : C —> C be a real entire function of exponential type at most 2 such that 
G € L 1 (R, |^(x)|- 2 dx). Then there exist W,T £ TL{E V ) such that 

G(z) = W(z)W*(z) - T{z)T*{z) 

for all z £ C. 

Proof. This is im Theorem 23]. □ 


5.2. Proof of Theorem O 
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5.2.1. Proof of parts (ii) and (iii). We shall prove these two parts together. Throughout this proof we denote 
by B{r) the closed ball of radius r centered at the origin, and the parameter v > — 1 is regarded as fixed. 

Existence. 

We prove first the existence of minorants and majorants that interpolate our given function at the ’’correct” 
radii. This is done in three steps. First we consider the particular dimension N v = [2i/+2]. A symmetrization 
argument will prove the existence for dimension N = 1 and, finally, the radial extension will allow us to 
conclude the result for arbitrary N. The reason to start with the dimension N v lies in the fact that, in this 
dimension, all considered functions are known to define tempered distributions, which allows us to employ 
the Paley-Wiener theorem for distributions to establish the exponential type of the extremal functions. Since 
symmetrization and radial extension preserve the exponential type, their application allows us to establish 
the required exponential type first for N = 1, and finally for arbitrary N. 

We consider below the minorant problem (the majorant is analogous). 

Step 1. The case N v = \2u + 2]. 

Recall that Q llt N u is assumed to be a tempered distribution satisfying 

/•OO 

G»,nA t)= / A-^/ 2 e-^l 2 /A d MA) (5-6) 

Jo 

for |t| > 5/2tt (this follows from ( 11 . 261 ) '). The function x D~ N (S,p,x) defined in ( 15 . 21 ) belongs to 
L 1 (M Ar,/ , |x| 2l " +2_Ar,y dx). In particular, this function is integrable in B( 1) (here it is important to have 
N v > 2v + 2) and thus defines a tempered distribution. 


Lemma 15. Let ip be a Schwartz function that vanishes in R(^- + e), where e > 0. Then 



„(<5,A,x) £(x) dx = 0. 


(5.7) 


We give the proof of Lemma [15] below. Assuming its validity, we show next that 

D~ Nu (6, p,t)= / A-^/ 2 e—l*l 2 / A d/r(A) 

for 111 > 5/2ir. In fact, if ip is a Schwartz function that vanishes in B{4^ + e), where e > 0, then 
[ D v n„ (^> Mi x ) ?( x ) dx = / [ |e _ '"' A|x|2 — (5, A, x)| </?(x) d/r(A) dx 

= [ [ |e“’ rA|x|2 - A,x)W(x)dxd/x(A) 

Jo Jwl n « k } 

= [ f X~ Nv/2 e~ n ^ 2/x ip(t) dt d/x(A) 

Jo Jr n v 

= f ( /°°A- Jv ^ 2 e-" |t| 2 /A d/r(A)^(t)dt. 

J\J 0 / 


(5.8) 


(5.9) 


From (15.61) and (15.81) we see that the function 

x Qn, Nv (x) - D~ Nu (d, n , x) 

defines a tempered distribution, with Fourier transform supported on B(5/2tt). It follows by the converse 
of the Paley-Wiener theorem for distributions that this function is equal almost everywhere on R' v " to the 
restriction of an entire function of exponential type <5. We call this entire function z i —> (<5, /z, z). By 
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Lemma [T3l and the assumption that x i —> is continuous on S. N ‘' \ {0}, we find that 

(6, n, x) = Nv (x) - D~ Nv ( S , p, x) 

for all x ^ 0. Since x D~ N (6, /x, x) is nonnegative, this implies that 

£v,n u (8, x ) < Gh,n v (x) 

for all x £ R ;V ", and by Lemma [13] it follows that 

{8, /x, x) = g,i,N v (x) 

whenever ^4^(<5|x|/2) = 0. An application of Fubini’s theorem gives us 

/ {£ M dV„( x ) - AvW.,(<*i Ab x )} | x | 2l '+ 2_JVl ' dx = f U^~(8,\)dp(X). 

J R N f Jo 

Hence, to conclude the proof of Theorem [3] for N = tV„ it remains to show Lemma fl5l 


(5.10) 

(5.11) 


(5.12) 


(5.13) 


Proof of Lemma\T5[ Recall from Theorem [5] that 

U v ,n v (5, A, z) = if Nt/ (L(Al, A, -))(8z/2), 


where z L(A 2 , A, z) is even, of exponential type at most 2, and belongs to L 1 (R, |x | 2l ' +1 da;). This implies 
that ^ i—^ L(A 2 ,A, z) must have at least one zero, for otherwise we would have L(A 2 1 X,z) = e cz+d with 
|c| < 2, contradicting the fact that it belongs to lA(R, |x | 2l ' +1 dx). Then, if a is a zero of 24 L(A 2 , A, z), so 
is —a since the function is even. Hence, J defined by 


J{z) = 


L(A 2 v ,X,z) 


(z 2 - a 2 ) 

is entire and has exponential type at most 2. It follows that 


J(z) := tf Nv (J){S z/2) 


{8 1 A, z) 

OT(zf + ... + 4)-a 2 


is entire of exponential type at most S. Since x i —> C v ,n v (■S , A, x) belongs to L 1 { , |x | 2l/+2 Nv dx), we have 
that x J{x) £ (here it is important to have 7V„ < 2v + 4). By a classical result of Plancherel 

and Polya [30], J is bounded on R^" and thus belongs to L 2 (R JV "). We now use the Paley-Wiener theorem 
in several variables (44] Chapter III, Theorem 4.9] to conclude that J is a continuous function supported on 
B(S/2tt) and hence, by Fourier inversion, 


J( z) = [_ J(t) e 2mt z clt. 

J B(5/ 2tt) 

This shows that |j7(z)| < Ce 5 ^, for z = x + iy, and thus 

(J, A, z)| <C(|<Iz/2| 2 -a 2 )e <5 l y l. 


(5.14) 


The bound (|5.14D allows us to invoke the Paley-Wiener theorem for distributions [28l Theorem 1.7.7] to 
conclude that the distributional Fourier transform of x h> (8, A, x) is supported on B(5/2n) and thus 
(15.71) holds. □ 


Step 2. The case N = 1. 

By (15.101) . the entire function z i —> {8, p, z) is radial when restricted to R^". This implies that 

(*5, /X, z) — ( 8 , p,z) 
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and Lemma [12] (ii) implies that z £ v ,n„ (8, p, z) admits a power series representation as in (14.41) . Define 
the entire function z £„(($, p, z) by 


L v {8, p, z) = £„ tN „(8,V,(z,0,...,0)). 


(5.15) 


By Lemma [T^l the function 2 K > L„(5, p, z) has exponential type at most <5, and we plainly observe that the 
analogues of (15.111) and (15.121) for N = 1 hold for all real x. A change to radial variables in (15.131) gives 


[ {gy.{x)-L v {8,p,x)}\x\ 2v+x dx = [ Ul (8, A) dp(\). 
Jr Jo 


(5.16) 


Step 3. The case of general N. 

We now define for z € C N 

£u,n( 8 , ix, z) = %j} N {L v {8 , p, -))(z) 

with L u (8,p,-) from (15.151) . By Lemma 1121 this function has exponential type at most 8, and we observe 
again that the analogues of (15.111) and (15.121) continue to hold. It follows from (15.161) and Theorem [2] that 


/ {S^a^x) - C v>N (S,p,x)} 

Jr n 


\2v+2 —JV dx _ 


pOO 

/ U^~(8, A) d/x(A). (5.17) 

Jo 


Optimality. 

Let C € (G^.n) be such that 

/ {£ m ,jv(x) - £(x)} |x| 2l/+2_Ar dx < oo. 

Jr n 


(5.18) 


By changing £ by £ we may assume that £ is radial on (note that this does not affect the value on the 
integral (15.181) 1. In this case, x >->• {C VjN (S, p, x) — £(x)} € £ 1 (R Ar , |x| 2 " +2-Ar clx). Let L : C —> C be defined 
by L{z) = C({z, 0,..., 0)). We want to show that 

/ {L v (8,fj,,x) - L(x)} \x\ 2 " +1 da; = f {£„,jv(6,/z,x) - £(x)} |x| 2l/+2_Ar dx > 0. (5.19) 

J — oo Jr n 

Let 

G(z) = L v {8, p, 2 z/8) — L(2z/8). 

Then G has exponential type at most 2, and by m it belongs to L 1 (R, \E V (x)| 2 dx). We are now able 
to use Lemma IT4l to write G = WW* — TT*, with W, T G 'H(E V ). It follows from (14.21) and the fact that 
{z K u (£, z); A v (£) = 0} is a complete orthogonal set in T-L(E U ) that 

/ OO 

{\W{x)\ 2 -\T(x)\ 2 }\x\ 2v+1 ds 

-oo 

poo 

[\W{x)\ 2 - \T(x)\ 2 } |.E„(aO|~ 2 dx 


= c„ 


= c v % Y 


\W(0\ 2 -\T(0\ 2 


A, (.0=0 




= c v 1 Y 


A„(0 =0 


g (0 


> 0 . 


In the last step we have used the interpolating property (|5.121) . which translates to 

L„(8, p, 2 t/S) = g„(M/8) > L(2t/8), 
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whenever A v (£) = 0. This establishes (15.191) . 


5.2.2. Proof of part (i). For n > 0, recall the definition of the measure given in (11.271) . We observe from 
part (ii) of Theorem [3] and from part (i) of Theorem [2] that 

pOO pOO poo 

U?~(6,ijl)= U?~(5,\)dn(\)= k 2v+2 U?-(k5,k 2 A) d/i(A) = k 2 "+ 2 / U?~(k5, A)d/z K - 2 (A) 

Jo Jo Jo 

= k 2i ' +2 U?-(k6,» k -2), 
and analogously for the majorant case. 

6. Applications 

6.1. Some examples. We start with a collection of interesting examples of radial functions included in our 
Gaussian subordination framework. 


6.1.1. Functions with exponential subordination. In d, the extremal problem (to - (El was solved for a 
class of functions of the form 

pOO 

G(x)=/ {e-^-e- T }do{r), (6.1) 

Jo 

where a is a nonnegative Borel measure on (0, oo) satisfying 


_j_ -y-2z/ + 3 


dcr(T) < oo 


( 6 . 2 ) 


for the minorant problem, and 


1 + T 


d(j(T ) < oo 


for the majorant problem. This includes examples such asxi->- log |x| and x i —> T(a + 1){|x| “ 1 — 1} for 
—2 < a < 2v + 1. Let us show that these functions are included in our Gaussian subordination framework. 

Let N = \2v + 2]. We start with the well-known identity 


d -i*i = 


y/u 


_M 1 , 
e 4 “ du 


and a change of variables gives us 


o - t | x | = 2_ f°° e -ir\\x\ 2 e -T 2 /(4n\) dX 

2 tt Jo 


A 3 / 2 ' 


Taking Fourier transforms on both sides leads to 

T T 

C N 


(t 2 + 47t 2 111 2 ) 2 2tt 


\~N/2 e -7r|t| 2 /A e -r 2 /(4^A) 


dA 
A 3 / 2 ’ 


where Cn = 2 N tt( n (see [H p. 6]). A standard argument as in (15.91) gives us 


[ G(x) v?(x) dx = [ \Cn f 
Jr n Jr n y Jo 


— - , 91 , |?N W + 1 d(T ( T ) ¥>(0 dt 

(T 2 + 47T 2 |t| 2 )— J 


X -N/2 e -n\t\ 2 /X 


R w \J 0 


0 -t 2 /(4ttA)_ 


27T 


dcr(T) 


dA 

A 3 / 2 


<p(t) dt 


for any Schwartz function <p : 
condition (Q4) with 


C vanishing in an open neighborhood of the origin. This establishes 

dA 


MX)= {1 e- 2 /( 4 ^da(r)) 


A 3 / 2 ' 


(6.3) 
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For the minor ant problem we must verify that, under (16.21) . the measure d/i defined by (16.31) verifies (11.231) . 
The proof for the majorant case will be analogous. In fact, 


pi r i 

Xdp(X) < 



dX 

Ab/2 < 00 ’ 


and 



X^+i MX) 




e ~T 2 /(4nX) 



dX 

^u+5/2 




-r 2 /(4n X) dA A 
\v+5/2 J 



< 




dA 




e -l/(4ny) 


dy \ 

yU+5/2 J 


T 


(27t) t 2i/ + 3 


dcr(r) 


< oo. 


Note that the Gaussian nx2 is not included in the class (16.11) . In fact, suppose that 



for x > 0, where the measure a satisfies (16.21) . Since 


'da 


< t e 


U — II 


(6.5) 


J i 

for Re (z) > 0, where c = min{Re (z), 1}, we see that both sides of (16.41) have analytic continuations to the 
half-plane Re (z) > 0. However, because of (16.51) . the right-hand side grows linearly when z = 1 + iy, y —> oo, 
while the left-hand side is given by e -7r ( 1+2ly-y ). 


The previous discussion shows that the class given by the Gaussian subordination method is strictly larger 
than the class (16.11) obtained in mi, showing that this becomes indeed the most powerful framework to deal 
with this sort of extremal problems. 


6.1.2. Positive definite functions. We might consider a finite measure /i and the function given by 

pOO 

9 n(x)= e -wAx 2 d/x(A). ( 6 . 6 ) 

Jo 

It is straightforward to check that this function fits our Gaussian subordination framework and thus we can 
solve the extremal problem for its radial extension to M. N , for any N > 1. Recall that, by a classical result 
of Schoenberg m Theorems 2 and 3], a function g admits the representation (16.61) if and only if its radial 
extension to 1 N is positive definite for all N > 1 or, equivalently, if the function x g(\/x ), for x > 0, is 
completely monotone (a function / : [ 0 ,oo) — > R is said to be completely monotone if (— 1 )” > 0 for 

t > 0 and n = 0,1, 2,..., and /(0) = /(0 + )). Among such functions we highlight: 

(i) x i—^ e - “l x l , for 0 < r < 2 ; 

(ii) x i-)- (|x | 2 + a 2 ) -/3 , for a, /3 > 0 ; 

(iii) x4 — log[(|x | 2 + a 2 )/(|x | 2 + j3 2 )] for 0 < a < fi. 

In fact, the Gaussian subordination method also works if we take a = 0 in the family (iii) above. The 
one-dimensional extremal problem for the family (iii) was particularly important in bounding the modulus 
of the Riemann zeta-function on the critical strip under the Riemann hypothesis [ICES]. 

6.1.3. Power functions. Let y(s) = 7r -s / 2 T(|). The function x i —> y(— a) |x| CT for a > —2v — 2 (minorant) 
and a > 0 (majorant), with a 0,2,4,..., also falls under the scope of our Gaussian subordination method. 
In fact, in dimension N v = \2v + 2], the distributional Fourier transform of this function outside the origin 
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is given by y(N v + a) |t 


—N u —(t 


(see [331 Chapter III, Theorem 5, Eq. (33)]) and we have 


/»oo 

7 (jV„ + <r) e ~ 7 r l t l 2 / A d/r CT (A) 

Jo 


(6.7) 


for dp a (X) = A 2 1 dA. It is clear that this measure satisfies (11.2311 (minorant) and (11.241) (majorant). 

6.2. Hilbert-type inequalities. The solution of our extremal problem is connected with inequalities of 
Hilbert-type. This was previously observed, for instance, in mmmm- For an account on the original 
Hilbert’s inequality and its weighted version, see the work of Montgomery and Vaughan [55] . In our setting, 
we let N > 1 be a given dimension, and consider the case 2^ + 2 = TV. For a nonnegative Borel measure p 
on ( 0 ,oo) we define 

/»00 

Q m ( y)= / X- N/2 e~^ 2 / x d^X). ( 6 . 8 ) 

Jo 

Note that the integrand is the Fourier transform of the Gaussian .Fa(x) = e —7 rA l x l in dimension N. 

Theorem 16. Let S > 0, N > 1 and 2v + 2 = N. Let {a,j}fLi be a sequence of complex numbers and 
{yj}jLi be a sequence of well-spaced vectors in R^, in the sense that |yj — yz| > <5 for any j ^ l. The 
following propositions hold: 

(i) If (i is a nonnegative Borel measure on (0, 00 ) satisfying (11.231) and Q^ is defined as in (16.81) then 


M M 

- u*-(2nd,p) Y \ a j\ 2 < Y j - y0- 

3 =1 


(6.9) 


3,1 = 1 


(ii) If n is a nonnegative Borel measure on (0, oo) satisfying (|1.24D and Q M is defined as in (| 6 . 8 |) then 


M M 

Y a o 2^(y j - yi) < u^ + (2nS, p) Y Kl 

3 =i 


( 6 . 10 ) 


3.1 = 1 


Proof. The minorant x i— > 27r<5, A, x) produced in Theorem [5] is integrable and of exponential type 27 t<5. 

Thus its Fourier transform is supported on B(8). We now observe that 

2 

,-27riy r x 


E a i e 

j =i 


/ OO 

{J'a(x) - £ V) jv(27r5,A,x)} 

-OO 

m />oo 

= Y a .i I {^a(x) - £„ tN (2nd, A,x)} e _2,rl(yj_yi )' x dx. 

• / i J —OO 


3,1 =1 


This leads us to 


M M 

- U?-{2TT6,X)Y\ a j\ 2 < Y a 3^iB(yj-yi)- 

3=1 


( 6 . 11 ) 


3,1 = 1 


Integrating both sides of (16.111) against dp(X), we arrive at (16.91) . The proof for the upper bound (16.101) is 
analogous, now using the majorant produced in Theorem [2] □ 

An interesting example arises when we consider the measures d/r CT (A) = A _ 2 _1 dA from 96.1.31 The 
identity (16.71) remains valid as long as cr > —N and Theorem [TO] gives us 

M M — 7 T jV+/V) c \ M 

a 3 ai ^ Uf + (2 tt 5, p a ) 


Uff (27 x5,p a ) 


n-{N+v)/2 r (A±Z) 


EM 2 i 


< 


El 


^ |yj - y i\ N+a - tt-(^+-)/ 2 r (^) ^ 


( 6 . 12 ) 
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The lower bound holds as long as a > —N, while the upper bound holds as long as a > 0. This inequality 
is related to the discrete Hardy-Littlewood-Sobolev inequality [24) p. 288], and improves the result of fill , 
where the restriction — N < cr < 1 had to be imposed. The one-dimensional version of (16.121) appeared in 
jTSl Corollary 22]. 

6.3. Periodic analogues. The interpolation tools developed in Section [2] allow us to treat the problem 
of optimal one-sided approximation of even periodic functions by trigonometric polynomials. This involves 
the same circle of ideas discussed in tm Section 8], in which the exponential replaces the Gaussian as the 
base function, and we refer to this work whenever possible. As in the Euclidean setting, the Gaussian 
subordination method also produces a strictly larger class in the periodic setting, when compared to the 
exponential method of m- 

We denote by ID C C the open unit disc and by <9D the unit circle. As in m, we identify measures on 
R/Z with measures on the unit circle c® via the map x i-» e 2nzx . Throughout this section we let 0 be a 
nontrivial even probability measure on R/Z (recall that a measure on R/Z is trivial if it has finite support), 
and thus the corresponding nontrivial probability measure on the unit circle c® (that we keep calling 9) 
satisfies 

0(A) = '0(A) 

for any Borel set A C c®, where A = (z : z £ A}. For Im (r) > 0, we consider the Jacobi theta-function 

OO 

9 3 (z,t)= J2 e RiTn2 +27rinz . 

n =—oo 

Poisson summation implies, for A > 0, that 

CO 

X-h^z,^ 1 ) = (6.13) 

j=—oo 

This is the periodization of the Gaussian T\[x) = e~ vXx2 . We first address the problem of majorizing 
and minorizing this function by trigonometric polynomials of a given degree n, in a way to minimize the 
id (R/Z, d$)-error. This problem with respect to the Lebesgue measure was treated in [12]. Here, of course, 
a trigonometric polynomial m(z) of degree at most n is a 1-periodic function of the form 

n 

m(z) = J2 a k ^ tikt , 

k——n 

where a*, € C. We say the m(z) is a real trigonometric polynomial if it is real for real z. Let us denote the 
space of trigonometric polynomials of degree at most n by A n . 

6.3.1. One-sided approximations to theta-functions. We first state the analogue of im Lemma 28]. 

Lemma 17. Let G be an even, 1-periodic, non-constant Laguerre-Polya function of exponential type t(G). 
Let z i ^ L{G, A, z) and z >->• M(G, A, z ) be defined as in Lemma 0 

(i) If G( 0) > 0, the function 

£(G,X,z) = J2l(G,X,z + j) 
jez 

is a real trigonometric polynomial in z of degree less than t(G) /(2tt) satisfying 

G(x)^X~^9 3 (^x,iX~ 1 ) — f(G,A,x)j > 0 


26 


(6.14) 




for all i€l, and 




for all £ £ R with G(£) = 0. 

(ii) If G has a double zero at the origin and G(0 + ) > 0, the function 

m(G, A, z) = yM(G, A ,z + j ) 
jez 

is a real trigonometric polynomial in z of degree less than t(G)/(2-7t) satisfying 

G(a;)|m(G, A, a;) — A _5 0 3 (x, *A^)| > 0 (6.15) 

for all i£l, and 

m(G,X^) = \-h 3 (^i\- 1 ) 

for all f £ R with G(£) = 0. 


Proof. The proof of m Lemma 28] uses that L has exponential type and is in L^R), that it minorizes 
x i—^ e -A l a: l on the real line, and that it interpolates at the zeros of G. When replacing the exponential with 
the Gaussian, these assumptions are given for L in our setting, hence the proof goes through without any 
changes; likewise for the majorant. □ 


In what follows, if an algebraic polynomial P(z) has degree n, we write P*{z) = z n P(z _1 ). Let {ip n }%L 0 
be the sequence of orthonormal polynomials in the unit circle with respect to the inner product 


{f, 9)L\aoA^)) = f{z)g(z)dd(z). 

Job 

We note that, since the measure D is even, the algebraic polynomials ip n have real coefficients (see [TT1 
Lemma 25]). Define the algebraic polynomials 

A n+ i{z) := ^*+1(2) + <p„+i (2)}, (6.16) 

Bn+l(z) := \{p*n+l (2) - <Pn+l(2)}, (6.17) 

which have degree n + 1 and have only simple zeros lying on the unit circle. Also, define the function 


K n (w,z ) = 


Pn+l (Z) Pn+1 O) - Pn+ 1 (2) (fin+1 (w) 


(6.18) 


1 — wz 

which is the reproducing kernel of the space V n of algebraic polynomials with degree at most n, with inner 
product (•, •)L 2 ( aD d ^) (for details see (TTJ Sections 8.1 and 8.2]). With these definitions we can state our next 
result. 


Theorem 18. Letn £ Z + andd be a nontrivial even probability measure on M/Z. Let ip n +i(z) = f n +1 (z; d$) 
be the (n + 1 )-th orthonormal polynomial in the unit circle with respect to this measure and consider 
K n , A n+ i, B n+ i as defined in (16.161) . (16.171) and (16.181) . Let A n +\ = {£ £ R/Z : A n+1 (e 2m ^) = 0} and 
B n+ 1 = £ R/Z : B n+1 (e 2 ™«) = 0}. 


(i) If t £ A„ satisfies 
for all x £ R then 


l{x) < A i0 3 (x,iX x ) 


£{x) d'd(x) < 


E 

£6-4n + l 


AT„(e 27ri «,e 27ri «)' 


27 


(6.19) 


( 6 . 20 ) 














Moreover, there exists a unique real trigonometric polynomial z >->• £#(n,X,z) £ A„ satisfying (16.191) 
for which the equality in (16.201) holds. 


(ii) If m £ A n satisfies 
for all x £ R then 


m(x) > A ^6^{x,iX 1 ) 

m{x) dd(x) > .„ —r-—. 

, w w ~ ^ K n (e 2 ™t,e 2 ™t) 


( 6 . 21 ) 


( 6 . 22 ) 


Moreover, there exists a unique real trigonometric polynomial z m$(n, X, z) £ A n satisfying (16.211) 
for which the equality in (16.221) holds. 


Proof. As is the case with the previous lemma, the proof of Theorem [TS] carries over verbatim from jlll 
Theorem 27], and we omit the details. □ 


6.3.2. One-sided approximations to general even periodic functions. Our goal now is to integrate the param¬ 
eter A against suitable measures and solve the external problem for a class of even periodic functions. We 
consider two classes of nonnegative Borel measures c on (0, 00 ). For the minorant problem we require that 

A _ 1 e _ wA - 1 -aA d <-( A ) < oo ( 6 _23) 

for any a > 0, whereas for the majorant problem we require the more restrictive condition 

(l + A - ^) e _7I ' A_1 dc(A) < 00 . (6.24) 

In this subsection we majorize and minorize the periodic function 

roo 

K(x) := J |A”^0 3 (x, iA -1 ) — A“50 3 (i ; jA -1 ) j dc(A) (6.25) 

by trigonometric polynomials of a given degree n, minimizing the L 1 (R/Z, d$)-error. Note the convenient 
subtraction of the term A“50 3 (I ; iA _1 ) to generate a better decay rate in A as A —> 0. If c satisfies (16.231) . h ? 
is well-defined and finite for all x Z (we may have h^( 0) = 00 ), and if satisfies (16.241) . h^ is well-defined 
and finite for all i£l. We obtain the following result. 




Theorem 19. Letn £ Z + andd be a nontrivial even probability measure onK/Z. Let ip n+ i(z) = <p n +1 
be the (n + 1 )-th orthonormal polynomial in the unit circle with respect to this measure and consider 
K n ,A n+ i,B n+ i as defined in (16.161) . (16.171) and (16.181) . Let A n +± = {£ £ M/Z : A n+1 (e 2nl ^) = 0} and 
Bn+1 = {£ e R/Z : B„+i(e 2 -«) = 0}. 


(i) Let q satisfy (16.231) . If £ £ A„ satisfies 


£{x) < h^(x) 


(6.26) 


for all x £ R then 

L- X, jM eilW (6 - 27) 

Moreover, there exists a unique real trigonometric polynomial z 1 —> £^(n,q,z) £ A n satisfying (16.261) 
for which the equality in (16.271) holds. 

(ii) Let c satisfy (16.241) . If m £ A n satisfies 


m{x) > h<;(x) 


28 


(6.28) 






























for all x £ R then 


/ m(x ) d$(:r) > 

ii/z 


E 

«eB„+i 


MO 

K n (e 2 ™Z, e 27ri «)' 


(6.29) 


Moreover, there exists a unique real trigonometric polynomial z i—> m$(n,q, z) £ A n satisfying (16.281) 
for which the equality in (|6.29l) holds. 


Proof. Part (i). Write 

h\(x) = X~^9 3 (x,iX~ 1 ) - X~^9 3 (^,iX~ 1 ). 

From Theorem 1181 the optimal trigonometric polynomial of degree at most n that minorizes h\ is 

U(n,X,z) := U(n,X,z) - X~^9 3 {\,iX~ 1 ). 


Let us write 

n 

U(n,X,z)= J2 Mn,X)e 2nikz , 

k——n 

where a k = a k ^. From the interpolation properties we have, for all A > 0, 

n 

7o(n, A, 0 = ]T o fc (n, A) e 2 " fc « = MO, (6.30) 

k=—n 

n 

U'{n, A,£) = E 27rifc a k (n, A) e 27rife ^ = h' x (f) (6.31) 

k——n 

for all £ £ A n +i (note that £ ^ 0). Since we have 2n+l coefficients {a k \ —n < k < n} and 2n + 2 equations, 
this is an overdetermined system. We also know that this interpolation problem has a unique solution, so 
we can drop the last equation and invert the coefficient matrix to obtain each a k as a function of the values 
{MO, MO; £ e M+i }. From the definition of h\ and identity (16.131) . together with (16.231) . we find that 


POO 

/ IMOIMA) 
J o 


< oo 


and 


IMOIMA) < oo 


(6.32) 


(6.33) 


for £ ^ 0. We then conclude that each coefficient a k (n , A) is absolutely integrable with respect to do(A) and 
we are able to define 

p OO 

X, z) dq(X). 


It is clear that 


0(n,?,x) < h q (x) 


for all x £ R and that 

= 0(0 

for all £ £ A n +i- The optimality and uniqueness of this minorant follows as in m Section 8]. 

Part (ii). Observe now that 0 £ B n + 1 and, when solving the system of equations for a k {n, A), we may throw 
away the equation corresponding to £ = 0 in (16.311) (which might not hold if h\ is not differentiable at the 
origin). Relation (16.321) continues to hold for all £ £ K./Z (due to the more restrictive condition (16.241) ') while 
(16.331) holds for £ ^ 0. The rest is analogous. □ 


Remark: In m Theorem 29] the analogous version of Theorem [T9] was proved for even periodic functions 
with an exponential subordination. Theorem [T9] is strictly more powerful than m Theorem 29] in the 
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sense that every function that fits the exponential subordination framework of El Theorem 29] also fits the 
Gaussian subordination framework of Theorem |T9l and there are functions (e.g. the theta-functions (16.131) 1 
that fit the Gaussian subordination framework but do not fit the exponential subordination framework of 
El Theorem 29]. Verifications of these claims follow the same circle of ideas described in 36.1.11 and are left 
to the interested reader. 
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